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Abstract 

Starting from a master equation in a quantum Hamiltonian form and a cou- 
pling to a heat bath we derive an evolution equation for a collective hopping 
process under the influence of a stochastic energy landscape. There results 
different equations in case of an arbitrary occupation number per lattice site 
or in a system under exclusion. Based on scaling arguments it will be demon- 
strated that both systems belong below the critical dimension d c to the same 
universality class leading to anomalous diffusion in the long time limit. The 
dynamical exponent z can be calculated by an e = d c — d expansion. Above 
the critical dimension we discuss the differences in the diffusion constant for 
sufficient high temperatures. For a random potential we find a higher mobil- 
ity for systems with exclusion. 
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I. INTRODUCTION 



Many systems behave on the phenomenological level essentially randomly and therefore other 
approaches for the theoretical treatment have to be employed. The randomness, resulting 
from stochastic forces or be intrinsic in the underlying microscopic theory, inevitably leads 
to the description of such systems in terms of probabilities and expectation values!'!. The 
time development of probability is usually found using a master equation. The past years 
have seen an exciting new development based on the observation! of the close relationship 
between the Markov generator of the master equation and a time evolution operator acting 
on a many-particle Fock space!'!, for some recent reviews compare!'!. The new insight has 
led to a series of remarkable exact solutions for the stochastic dynamics of interacting particle 
systems, for a recent overview see!. Despite of exact results the mentioned method has been 
also fruitful in an approximative description of other models such as the facilitated kinetic 
Ising system as a candidate for glassy systems!~0 or in branching and annihilation random 
wa lks@. Whereas the original paper!, see alscM are concerned with a mapping of the master 
equation to a representation in terms of second-quantized bosonic operators a great progress 
for exact solvable models had been achieved by mapping to spin-1/2 Pauli-operatorJl. This 
mapping to spin systems applies to processes where each lattice site can be occupied by only 
a finite number of particles. Physically, this restriction may be hard-core constraints or fast 
on-site annihilation processes. Obviously, such a mapping simulates the exclusion principle 
for classical lattice models with in a cellular automata. 

In the present paper the Fock space description is applied for systems far from equilibrium 
which are coupled to a heat bath. In particular, we discuss the collective hopping process of a 
classical many body system, coupled to the mentioned heat bath, and under the influence of 
a random energy landscape realized by a stochastic activation energy. The particles making 
random walks have to overcome spatially distributed energy barriers. As the consequence the 
hopping process is accomplished by a competing force field which can give rise to anomalous 
diffusion. Further, the analysis should be different considering both cases, the bosonic and 
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the exclusive ones. In the first case the particles should find more rapidly the local energy 
minima however because of that their mobility could be reduced. As the consequence of 
the random walk where the particles have to overcome spatially distributed energy barriers, 
the resulting effective force field can give rise to anomalous diffusion. It is well known 
that one of the reasons for an anomalous diffusive behavior can be traced back to the 
influence of a stochastic force field below a critical dimensionBEl. An alternative way of self 
induced anomalous diffusion had been discussed recently^ introducing a feedback coupling 
between the diffusive particle and its local environment. Both, the disorder and the memory 
controlled feedback may lead to a subdiffusive behavior or to localization. The analytical 
approachEl could be confirmed by simulations in one and two dimensions where at the 
critical dimension d c = 2 logarithmic corrections in the mean square displacement had been 
foundSfi 

Here, we demonstrate that the Fock space approach leads in both cases, bosonic and under 
exclusion, may lead to anomalous diffusion. Within the long time limit and on a large spatial 
scale both systems belong to the same universality class. 



II. QUANTUM APPROACH TO NONEQUILIBRIUM 

The analysis is based on a master equation 

d t P(n,t)=L'P(n,t) (1) 

where P(n,t) is the probability that a certain configuration characterized by a state vector 
n = (jii, ri2 ■ ■ ■ tin) is realized at time t. There are two special cases, eithe each lattice 
site is occupied by an arbitrary number of particles Hi = 0, 1, 2 ... or as in a lattice gas 
rti = 0, 1. Further, the occupation numbers rij are considered as the eigenvalues of the 
particle number operator defined by creation operators d\ or by annihilation operators di. 
The problem is to formulate the dynamics in such a way that the possible realizations for 
the occupation numbers are taken into account explicitly. The situation in mind can be 
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analyzed in a seemingly compact form using the master equation in a quantum Hamilton 
formalismlili'El, for a recent reviews seeol. The dynamics is determined completely by the 
form of the evolution operator L', specified below, and the commutation relations of the 
underlying operators d\ and dj. Within that approach^ the probability distribution P(n,t) 
is related to a state vector | F(t)) in a Fock-space according to P(n,t) = {n | F(t)). The 
basic vectors | n) are composed of the operators d\ and d{. Using the relation 

\F(t)) = J2P(n,t)\n) (2) 

rii 

the master eq. ([!]) can be transformed into an equivalent one in a Fock-space 

d t | F(t)) = L | F(t)) (3) 

where the operator L' in ([!]) is mapped onto the operator L = J2 \ ™)L' mn {n | in eq.(||). 
It should be emphasized that the procedure is up to now independent on the realization 
of the basic vectors. Originally, the method had been applied for the Bose caseiQ0. 
Recently, an extension to restricted occupation numbers (two discrete orientations) was 
proposed&ilii. Further extensions to p-fold occupation numbers^ as well as to models with 
kinetic constraints^ and to systems with two heat bathesil are possible. 
As shown by Doii the average of an arbitrary physical quantity B(n) can be calculated by 
the average of the corresponding operator B(t) 

(B(t)) = J2P(n,t)B(n) = (s \ B \ F(t)) (4) 

rii 

with the state function (s |= Y^,{n I- The evolution equation for an operator B(t) reads now 

d t {B) = {s\[B(t),L]\F(t)) (5) 

As the result of the procedure, all the dynamical equations governed by the classical problem 
are determined by the structure of the evolution operator L and the commutation rules of 
the operators. 
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III. COUPLING TO A HEAT BATH 



The evolution operator for a collective hopping process is different for an arbitrary occupa- 
tion number, denoted as Bose case, or an restricted occupation number, denoted as Fermi 
case. For the last system the operator Lf reads! 

L f = vYl ( d l d i -i 1 - n i) n j) ( 6 ) 

where /i is the hopping rate between adjacent sites % and j. The occupation number operator 
rij = d\di is related to the spin operator by the relation Si = 1 — 2nj and the commutation 
rule is [di, d 3 ] = 5^(1 — Irii). 
For the Bose case we get 

L b = vJ2{ d i d j -%) ( 7 ) 

i,i 

where d\ and di fulfills the Bose commutation rules. A generalization to processes under the 
coupling to a heat bath with a fixed temperature T is discussed inl. As demonstrated inll 
the evolution operator has to be replaced by 

L f = fi [C 1 - d i d )) exp(-/3H/2)4dj exp(j3H/2)] (8) 

where the hopping rate \x defines a microscopic time scale; (3 = T _1 is the inverse temper- 
ature of the heat bath and H is the Hamiltonian as a measure for the energy. A further 
generalization is realized by introducing different local heat bathes is discussed in0. In the 
bosonic case the generalization to finite temperatures leads to 

L b = ^[(l- <y exp(-(3H/2)d\d, exp(j3H/2)] (9) 

Here we study the case that the Hamiltonian H in eqs.(|],[D is simply given by a stochastic 
energy landscape defined by the energy functional 

H = J2e l n t (10) 
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Whenever the energy is positive the empty site is energetically favored. Further, e is assumed 
to be a stochastic local energy the distribution of which will be introduced below based on the 
continuous representation. In this manner, the model describes a collective hopping process 
where the jumping particles are subjected to a local random energy £j which supports or 
prevents the hopping process with a probability proportional to exp(±£j/2T). Taking into 
account the commutation rules we get in both cases 

e-f^dtef^ = d^' 2T e^ H ' 2 ' d\^ H ' 2 = d\e~^ 2T (11) 

Using eq. fl5|) and the algebraic properties of Pauli-operators, the evolution equation for the 
averaged density reads 

^ l d t (n r ) = YH ex P i( £ j ~ £ r)/2T) (rij) - exp ((e r - £j)/2T) (n r ) 

-2sinh(^^)(n r n,) (12) 

In the Bose case the evolution equation is much simpler. 

fi-%(nr) = £ [exp(( £ ,- - s r )/2T)( nj ) - exp((e r - e,)/2T)(n r )} (13) 

Both equations reflect the conservation of the particle number which will be more trans- 
parent in a continuum representation. In the special case of a constant energy e r = Ej it 
results the conventional diffusion equation in a discrete version. When the energy changes 
from site to site the nonlinear eq.(|l2"l) is the first step of a whole hierarchy of evolution 
equations. Assuming now smoothly changing energy e r and density n r a gradient expansion 
is appropriate up to the order I 2 where 1 is the lattice size. To make the expansion invariant 
under the underlying rotational symmetry we have to use the following identity 

5>xp(( £j - e r )/2T)(n j ) = ^(n r ) 

j(r) j(r) 

+ exp(-e r /2T) b x P{£j/2T)( nj ) - exp{e r /2T)(n r )} (14) 
Such an expression reads in a continuous representation including terms of the order I 2 

zn{r, t) + exp(-e(r) /2T) V 2 [(exp(e(r) /2T)n(r, *)] 
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with the averaged density (n r ) = n(r,t); z is the number of nearest neighbors. After 
decoupling the nonlinear term in eq.(^) and performing the continuous limit the density 
n(r, t) obeys the following nonlinear diffusion- like equation 

fTH~ 2 dtn = V 2 n + n(l - n)— + (1 - 2n)Vn ■ Ve/T (15) 

In a system with exclusion the density couples in a nonlinear manner to the stochastic energy 
field e(r). Due to the exchange coupling of the evolution operator L in eq.@ the resulting 
equation fll5"D is a conserving one where the current is given by 

if = -Vn-n(l -n)^ (16) 

In the Bose case we find after performing the continuous limit the density n(r, t) obeys the 
following exact equation 

ir x r % d t n = V 2 n + ^V[nVe] (17) 

The conservation law is manifested in the current 

j 6 = -Vn-n^ (18) 

The resulting equation is nothing else as the conventional diffusion equation under an ad- 
ditional drift term where the Einstein relation is automatically fulfilled. Remark that one 
can derive a similar equation when the system is coupled to two heat bathes with differ- 
ent temperatures. In that case one has to replace e(r)/T by where v is the chemical 
potential and T(r) is the local temperature, see also!. In the Bose case eq.(|l7D depends 
on the density in a linear manner. It is of Fokker-Planck-type when the density n(r, t) is 
considered as the single probability distribution to find a particle at site r at time t. Such an 
interpretation is always possible because we have not taken into account any interactions. 
Therefore, the particles are independent from each other and the concentration field behaves 
as the probability distribution of a single particle of this system. Different to the case of an 
arbitrary occupation the current jj includes a term n(l — n) which is characteristic for sys- 
tems with exclusion. Due to the exclusion principle the systems reveals a kind of correlation 



which leads even in the mean field limit to a nonlinear current. Following the discussion 
for the Bose case eq.(|T5|) can be interpreted as a nonlinear Fokker-Planck-equation for a 
single particle. The nonlinearity reflects a feedback of a particle to itself due to the excluded 
volume effect. 

It seems to be more appropriate to introduce the force vector f(r) = — Ve(r) the evolution 
equation in the Bose case reads now 

^-H- 2 dn(r, t) = V 2 n - • Vn - -V • fn (19) 

In the Fermi case the corresponding equation is 

// _1 r 2 5n(r, t) = V 2 n - -f ■ Vn(l - 2n) - -n(l - n) V - f (20) 

When the force field f(r) is a stochastic one the system offers anomalous diffusive 
behaviorBEl. 



IV. SCALING 

Now let us discuss both equations when the force field is an stochastic pure spatial dependent 
field, the correlator of which is given by 

/ a (r)/(r') = <Mr-r'), W) = (21) 

After averaging over the distribution function of the force field the system is homogeneous 
depending only on the difference of the spatial coordinates. The most general form of the 
function <ft ai is given in a Fourier representation by 

0o, 7 = A(q)(8 ai — n a n 7 ) + 5(g)n a n 7 with n a = — (22) 

Introducing dimensionless variables x — > aA _1 , t — > tA~ z , where z is the dynamical critical 
exponent and further n — > nA d and according to eq.Q2"^) for constant A and B f — > fA d / 2 we 
find the critical dimensionality d c — 2. For d < 2 the term proportional to V(fn) is relevant 
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whereas the additional term in case of exclusive motion oc nfVn is only relevant for d < 2/3. 
That means for the physical dimension d > 1 both models belong to the same universality 
class, where only d < 2 the disorder is relevant. Physically the result is obvious because 
in the long time limit and for a large spatial scale the Fermi system can be considered to 
consist of blocks with an increasing size. The larger such a block the more irrelevant is to 
distinguish both cases, arbitrary occupation and restricted occupation. In case of d < 2 the 
system reveals anomalous diffusive behavior as it had been demonstrated for a similar model 
not for the density n(r, t) but for the probability to P to find a particle at time t at the point 
r. Making the same calculation we end up with the flow equations for the dimensionless 
coupling parameters D = fxl 2 , a = D ^ T2 K d , b = D 2 T2 K d , with K d {2-K) d : the volume of 



the d-dimensional unit sphere and e = 2 — d £ = ln(^- 



D 



3D 
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d 
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d 



a 

e ~d 



(23) 



In the same manner one can derive an equation for the mean square displacement R 
A 2 s(D, a, b) with s = (r 2 ) The flow equation can be written as 



(24) 



That equation leads to a scaling behavior of the mean square displacement in the vicinity 
of the fixed points of eqs.(|2"3"|). In order to keep the diffusivity D fixed to its bare value the 
effective dynamical exponent z(£) satisfies z(£) = 2 + b{^)/d + a(£)(l — d)/d. When the 
disorder is irrelevant the fixed points are a* = b* = the exponent is z = 2. For the fixed 
point a* = ed, b* = it results z = 2 — e and for a* = b* = ed we find z = 2 + 0(e 2 ). These 
values are well known0'@. At the critical dimension d c = 2 we proceed on the following 
manner. The observation time t is related to an initial time to by 



t = t exp{[ z(£')d£') 
Jo 



(25) 



Using eqs.(|I||4] we can fix the scaling parameter £ according to eq.(^ to be 

^ 2 V 2 v 2 t 



where ao is initial value for the parameter a. From eq.(|2l|) we find the following behavior 
for the mean square displacement 

(r 2 ) = Cl f + c 2 f ln(l) (26) 

where ci and c 2 are two non-universal constants. As expected the system reveals logarithmic 
corrections at the critical dimension. 

V. BEHAVIOR ABOVE THE CRITICAL DIMENSION 

The thermalized version of the Fock space representation, see eqs.(|8],|9]), leads in the limit 
T — > oo to conventional diffusion. In the high temperature limit the particles are able to 
overcome each barrier and as the consequence of the stochastic hopping process one finds 
diffusive behavior in the long time limit independently on the underlying statistics. When 
the temperature is finite there appears a competition between two processes resulting in a 
different behavior for both systems. Bose particles can easily find a minimum within the 
energy landscape defined by the stochastic force. Particles with exclusion have to search for 
a longer time and on a larger scale to reach an appropriate potential minimum. From here 
one would conclude to an enhanced diffusivity. On the other hand, the mobility of Bosons is 
eventually reduced because they find more rapid a stable minimum. Due to the established 
universality for low dimensions a variation of the behavior should be only observed above the 
critical dimension. In this regime conventional perturbation theory should be applicable. 
Let us therefore present lowest order corrections to the the diffusion parameter D. The 
effective diffusivity is defined by 

D <»-fW\ . . < 27 » 

^ q=0,uj=0 
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As well as the Bose and the Fermi system lead in second order, proportional to J^, to non 



trivial corrections which are also manifested in the averaged density n(r, t) or the averaged 



correlation function n(r, t)n(r', t'). Indeed, the Fermi system offers additional terms for the 
density or the correlation function compared with the Bose case. However those terms does 
not contribute at zero wave vector and hence there are not relevant corrections to the di- 
vergent part of D e ff for d < d c . Above d c the behavior of the effective diffusion coefficient 
can be estimated using a perturbative approach around the homogeneous solution denoted 
by n. We get 

U eff ~ U eff + L 

with I = — f-I x [B - A{d - 1)] (28) 

Ii can be expressed by a momentum integral which is always positive in the mesoscopic 
regime A > I. For B — A(d — 1) > 0, realized for a pure potential field (B is the relevant 
variable, see eq.(|22|)), eq.(|28|) leads to 

D lff > D e.ff ( 29 ) 

Remark that the correction to the bare diffusion coefficient D is of the order (1 — 2n) 2 , that 
means for the half-filled case there are no corrections. That reasonable result should be also 
valid in a more refined approach. 

Because the homogeneous solution is not necessary a stable one we can also estimate the 
behavior using linear stability analysis around the stationary solution denoted as n s (r). Let 
us introduce n(r,t) = n s (r) + y(r,t) then the correction y(r,t) fulfills in the Bose case the 
equation 

d t y = DV 2 y + jV{yVe b ) with f (r) = -Ve 6 (r) (30) 

Here £&( r ) = e(r)—v is the true stochastic potential introduced by eq.fllQl) and v plays the role 
of the chemical potential which regulates the occupation number. In case of the exclusion 
model the deviation from the stationary solution y(r, t) satisfies the same equation however 
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one has to replace the potential in the Bose case, given in eq. (|30|) , by another effective 
potential 

cosh ( £ ^r~" ) 



e 6 (r) ->e/(r) = 2Tln 



cosh ^ 



(31) 



We have gauged the potentials so that for £/(r) = also e(r) = 0. The hopping particles 
under exclusion are subjected to the modified stochastic energy landscape given by £/. 
Expanding £/ in terms of e we find the relation 

£/ (r)~-tanh(^) £ (r) (32) 

From here it results 



£/ (r)e / (0)^tanh 2 (^) s b (r)e b (0) (33) 

The effective correlator of the disorder in the Fermi system is drastically decreased in com- 
parison to the Bose case. This result is compatible with the previous discussion leading to 
eqs.(^,^). In particular in the vicinity of half-filling (where the chemical potential v is 
zero) the influence of the disorder is very weak. This special case corresponds to vanishing 
linear term expanding eq.([3~Tl) according to powers of e. In the leading order we obtain 



s/(r) 



e 2 (r) 



AT 



Different to the Bose case the effective stochastic potential Ej, eq. (|31[) , is always positive 
definite, that means all the deep negative minima of the original stochastic potential become 
maxima and therefore they are not more available in case of Fermi system. Obviously, they 
are already occupied and hence they are not accessible for particles. 



VI. CONCLUSIONS 

In the present paper the collective hopping process on a lattice is studied systematically 
when the particles are subjected to a random energetic landscape manifested by a stochastic 
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energy profile. In particular, we have taken into account both cases, each lattice site is only 
occupied by one particle or each site can absorb an arbitrary number of particles. Physically, 
one expects different behavior. Whereas in the situation under exclusion a particle should 
spend more time for searching an appropriate energy minimum within the stochastic energy 
the bosons tend to reduce their mobility because they remain for a longer time in the 
local minima. A further influence on the motion of the particles is given by the coupling 
to a heat bath which supports the tendency that the system equilibrates. Starting on a 
master equation in a second quantized form both cases can be easily realized in terms 
of Bose-operators or spin- 1/2 Pauli-operators. The annihilation and creation process of 
particles leads in both cases to a density gradient characteristic for a random walk. Due 
to the additional coupling to a stochastic energy each particle can not follow that gradient 
simply but it has to overcome an energy barrier at its starting point and at its end point. 
There appears a conflicting situation that a particle follows the density gradient but the 
energy at the starting point is higher than at the end point. In this manner it will jump 
from an occupied to an empty site however under mobilizing a higher amount of energy 
(lower temperature). The other situation consists of the fact that a particle follows the 
density gradient and the energy barrier at the starting point is lower than at the end point 
(high temperature regime). In this case the hopping process is highly supported by the 
energy landscape whereas in the previous one the process is restricted. As the consequence 
anomalous diffusive behavior should be realized below the critical dimension. 
In the paper we have demonstrated that the Bose-as well as the Fermi-system belong below 
the critical dimension to the same universality class within the long time limit and on a large 
spatial scale. For an increasing scale the system can be considered consisting of blocks with 
an increasing number of particles. Thus, the cases of restricted and unrestricted occupation 
number per lattice site should be irrelevant. Despite of the universality the density and the 
correlation function of both systems are different, in particular for an intermediate interval. 
In particular, we have discussed the situation above the critical dimension where the diffusion 
constant can offer different behavior in both cases. 
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